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NOMENCLATURE 


matrix for rotation about the y-axis 
matrix for rotation about the z-axls 
wing span, m 
chord, m 

mean aerodynamic chord, m 

propeller or nacelle incidence (body centerline reference) , deg 
propeller or nacelle yaw angle, deg (see fig. 1(b)) 

Mach number 

nacelle with simulated jet exhaust 
propeller 

radius of propeller disc, m 
radial distance along propeller blade, m 
airfoil thickness, m 
column vector (see appendix A) 
velocity in the x-direction, m/sec 
components of the U-column vector 
transformed column vector, V = AU 
free- st ream velocity, m/sec 
velocity in the y-directlon, m/sec 
(see fig. 4) 

components of the V-column vector 
transformed column vector, W = BV 
rectangular wing 
swept wing 

tapered wing with a crank trailing edge 


*< 4 1 
• 

i 

m 

* 

W 4 twisted and cambered wing 

v 

w velocity in the z-direction, m/sec 

Wi,W 2 ,W 3 components of the W-column vector 
x,y,z Cartesian coordinates 

a angle of attack, deg 

<x T propeller blade angle of attack at — » 0.75, deg 

L K 

IT * 

3 propeller blade pitch angle at ■=• * 0.75, deg 

A T difference between maximum and minimum values of a. 

aL L 

5 j / azimuth angle, deg (see fig. 3) 

ij) (see fig. 4) 

a) rotational velocity, rad/sec 

?,n,5 (see fig. 3) 
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INTERFERENCE EFFECTS OF AIRCRAFT COMPONENTS ON THE 
LOCAL F/LADE ANGLE OF ATTACK OF A WING-MOUNTED PROPELLER 

J. P. Mendoza 
Ames Research Center 


SUMMARY 


A brief theoretical study was conducted at M = 0.6 to obtain an under- 
standing of the aerodynamic interference effects on a propeller operating in 
the presence of different wing-body-nacelle combinations. Tbs study was 
directed toward minimizing the unsteady blade angle-of-attack variation with 
azimuth angle by varying the pitch and yaw of the nacelle. For the particular 
configuration of interest the minimum blade angle-of-attack variation occurred 
with the nacelle pitched downward 4.5° and yawed inward 3.0°. 


INTRODUCTION 


Since 1973, the fuel fraction of the direct operating cost for air trans- 
ports has been steadily increasing, thus creating the need for fuel-efficient 
airplanes (ref. 1). Studies have Indicated that a turboprop-powered airplane 
operating at M = 0.8 could achieve a 10-20% savings in fuel relative to a 
comparable turbofan-powered airplane. For this reason, research efforts are 
currently underway in categories such as advance propellers, propeller noise 
and fuselage noise attenuation, propeller and gearbox maintenance, and 
airframe-propulsion systems integration. In this last category, both theoret- 
ical (ref. 2) and experimental (ref. 3) investigations have been conducted to 
determine the propeller slipstream effects on wing-body-nacelle- and wing- 
body combinations, respectively. One aspect of airframe-propulsion systems 
Integration that has not Ran widely investigated is the problem of the inter- 
ference effects on the propeller blade attributed to the presence of airplane 
components such as wings and bodies. In particular, the problem that has not 
been previously addressed is that of minimizing the cyclic bending moments 
applied to the propeller blade caused by the local blade angle-of-attack 
variation with azimuth angle. As a result, the present investigation was 
conducted: (1) to obtain a better understanding of the interference effects 

on the propeller blade due to the presence of wings and bodies and (2) to 
minimize the blade angle-of-attack variation with azimuth angle for a given 
turboprop transport model. 


AIRPLANE COMPONENTS 


The five different configurations used in the present theoretical study 
are shown in figures 1 and 2, They include an isolated nacelle with a simu- 
lated jet exhaust (fig. 1 ) and four different wlng-body-nacelle configurations, 
also with simulated jet exhausts (figs. 2(a) through 2(c)), As noted in 
figure 2(c), two of the configurations (PBW 3 N and PBW 4 N) were identical except 
for wing camber and twist. The wing sections for the configurations with the 
rectangular and swept wings had the same thickness distribution. The airfoil 
coordinates are presented in table 1, The coordinates for the tapered wing 
with the crank trailing edge are shown in table 2 for four span stations. 

The coordinates at four span stations for the cambered and twisted wing which 
had the same planform as the tapered wing are presented in table 3. Each of 
the four wings had 2° of dihedral. The nacelle was pitched and yawed about 
the fixed reference point shown in figure 1 . 


METHOD 


Because a generalized method is not presently available, an approximate 
method was developed for estimating the interference effects of nearby air- 
plane components on the local angle of attack of a propeller blade. The 
method is based on the assumption that the inflow into the propeller disc is 
dominated by the aircraft configuration and is essentially independent of the 
propeller and its slipstream. Under this assumption, the local Inflow veloc- 
ities can be combined vectorially with the rotational velocity of the propel- 
ler blade to define a local blade angle of attack as a function of azimuth 
angle. The method used to predict the local flow velocities was the Douglas- 
Neumann Potential Flow Program (ref- 4) which is a linear panel method capable 
of analyzing complete aircraft configurations. Using this method, velocities 
are computed at off-body points corresponding to points at r/R =0.75 on 
the propeller blade at different azimuth angles. The point at r/R = 0.75 
coincides with the centroid of the load distribution of the propeller blade 
and the flow at this point is considered to be representative of that for the 
entire blade. 

The problem of minimizing the cylic bending moments of the propeller 
blade caused by the variation in the local angle of attack of the blade is a 
difficult problem in itself. The difficulty is increased at higher subsonic 
Mach numbers where transonic effects are present and no adequate transonic 
analysis is presently available. To simplify the problem and allow the use of 
linear methods, it was assummed that the local angle-of-attack variation of 
the propeller blade at M = 0.8 is essentially the same as that at M = 0.6 
for the same velocity ratio which is the ratio of the tip velocity to the free- 
stream velocity. The velocity ratio was 1.0. The design blade pitch angle, 

8 0 75 R, of 56.5° at M = 0.6 was used throughout the present study. 

Shown in figures 3(a) and 3(b) are the flow velocities at a. point on the 
propeller blade. The propeller, unless otherwise noted, has right-hand 
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rotation (counterclockwise as seen by an observer In front of the airplane) 
and Is Installed on the right wing panel. The flow velocities which are 
computed by the method of reference 4 are transformed by the procedure 
described In appendix A Into the wg , and W 3 components shown In 
figure 3(a). These components, In turn, are resolved Into velocity components 
along the axes of a coordinate system that rotates with the propeller. Two of 
the three components are shown In figure 3(b). The third component that Is 
parallel to the radius of the propeller does not contribute to the blade bend- 
ing moment and, therefore, Is not Included In the analysis. From figure 3(b) 

v^ « rto - W 2 sin tj> - W 3 cos i|) 


Since 


<j> - tan“ 1 (w 1 /v N ) 

then c»l, the local angle of attack of the blade is given by 

« L - 3 - <f> 

where (3 is the propeller pitch angle. 


RESULTS 


Component Buildup 

To obtain a better understanding of the interference effects on the 
propeller blade attributed to the presence of nearby airplane components such 
as wings and bodies, an airplane component buildup was conducted starting with 
an Isolated propeller and continuing on to wlng-body-nacelle configurations 
with varying wing geometry. Blade angle-of-attack variations with azimuth 
angle were compared for the different configurations. By using the results of 
the Isolated propeller study as a basis for comparison, the effects of adding 
or changing various airplane components can be assessed. The local angle of 
attack of the propeller blade is understood to be computed at r/R = 0.75. 

Figure 4 shows the variation of the local angle of attack of the propeller 
blade with azimuth angle for an isolated propeller in a uniform flow field. 

The solid line represents the condition where the propeller axis of rotation 
is aligned with the free-stream velocity vector, while the dashed line repre- 
sents the condition where the propeller axis of rotation is pitched upward 2 ° 
which Is observed to produce a Acxl of 2.5°. 

Figure 5 shows the results for an isolated propeller P and for a pro- 
peller in the presence of a nacelle with a simulated jet exhaust PN. In both 
cases the propeller axis is at ia = 0° . The asymmetry of the nacelle induces 
nearly 1° of unsteady blade angle-of-attack variation. Figure 6 shows the 
effects of pitch angle on the PN configuration. Note the variation in 
AaL with varying i a . The smallest value of Aoil is at let = -0.5°. 
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Comparison of the blade angle-of-attack characteristics for the con- 
figuration buildup is shown in figure 7 starting with an Isolated nacelle and 
continuing on to a rectangular wlng-body-nacelle configurations PEW^N. For 
this comparison the body and wing are at 2° angle of attack and the nacelle 
is pitched downward 2.5° relative to the body centerline (-0.5° relative to 
the free stream). Note the small contribution of the body to the overall level 
of ai which is in sharp contrast to the effect due to the wing with its 
attendant upwash field. 

Figure 8 shows the effects of varying nacelle pitch angle on the blade 
angle-of-attack characteristics for the rectangular wlng-body-nacelle 
configuration, PBWiN. The wing and body are at 2° angle of attack while 
the nacelle is pitched from -2.5° to -4.0°. The smallest value of to j, 
occurs at 1$ * -3.5°. 

Figure 9 shows the effect of wing sweep. Blade angle-of-attack character- 
istics for wlng-body-nacelle configurations with a rectangular wing (PBWiN) 
and a swept wing (PBWaN) are compared. The wings and bodies are at 2° angle 
of attack and the nacelles are at ** -3.5°. The sweep angle for the swept 
wing was 35°. Wing sweep is shown to produce a substantial Increase in 
Aai because of the sidewash that is induced by a wing sweep. To compensate 
for the effects of sidewash induced by wing sweep, the nacelle for the swept- 
wlng configuration (PBW2N) was yawed from 0° to 2.5°, The results are shown 
in figure 10. The wing and body are at 2° angle of attack and the nacelle is 
at ig = — 3.5 P . The smallest value of Acxl is at ig * 2°. 

The effects of wing planform on the blade angle-of-attack characteristics 
were investigated using the swept wing-body-nacelle configurations PBWgN and 
PBW3N. The results that are shown in figure 11 show a small change in the 
blade angle-of-attack characteristics as a result of the change in wing plan- 
form. A comparison of the blade angle-of-attack characteristics for the 
tapered wlng-body-nacelle configuration with and without camber and twist is 
shown in figure 12. The significant changes shown in the blade angle-of-attack 
characteristics for the cambered and twisted wing are produced by the change 
in the Induced upwash field of the wing. 


Blade Angle-of-Attack Minimization 

In the present investigation, the procedure used to minimize the cyclic, 
bending moments applied to the propeller blades of a turboprop transport model 
is to minimize Aai. Except for the addition of nacelles and simulated jet 
exhausts, the PBW4N configuration is the same as that used in the investiga- 
tion reported in reference 3. Since it has been shown that A«l can h e 
minimized by varying the pitch and/or yaw of the nacelle, the nacelle of the 
PBW4N configuration was yawed from 2° to 3.5° in 0.5° increments. At each 
yaw angle the nacelle was pitched from -2.5° to -5.5° in 1° increments. The 
results of this study are shown in figures 13(a) through 13(d) and the data 
for these figures are summarized in figure 14. This shows values of Aa^ for 
each combination of pitch and yaw angles. The minimum value was found to be 
2° and corresponds to i a = -4.5° and ig = 3.0°. 
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To determine the effect of reverse propeller rotation which corresponds 
to a propeller with right-hand rotation mounted on the left wing panel, blade 
angle-of-attack characteristics for the PBW 4 N configuration are compared for 
the propellers with counterclockwise and clockwise (reverse) rotations. The 
body and wing are at ct * 2 ° and the nacelle is at l a * -4.5* and ig » 3.0° 
which are the optimum pitch and yaw angles for minimum Aaj, for the counter- 
clockwise rotating propeller, Figure 15 shows that in addition to the expected 
change in phase angle there is an Increase in Ackl from 2° to 3° (reverse 
rotation) . 


CONCLUSIONS 


The interference effects on the propeller attributed to the presence of 
different airplane components such as wings and bodies (including nacelles 
with/ simulated jet exhausts) were found to affect the blade angle-of-attack 
characteristics significantly. Compared to the effect of varying the inclina- 
tion of the propeller axis of rotation, however, these effects are not as 
large. Each component is shown to affect the blade angle of attack to some 
extent. The largest component effect came from the wing. The minimum value 
of ActL for the PBW 4 N configuration was obtained with a nacelle orientation 
of ! a “ -4.5° and ig = 3.0°. 


APPENDIX A 


As previously described, the nacelle can be pitched and yawed about a 
fixed reference point (fig, 1). For given values of % and lg, velocities 
can be computed (using the method of reference 4) at off-body points corre- 
sponding to points on the propeller blade, To compute the local blade angle 
of attack, these velocities are resolved into components along the axes of a 
rotating orthogonal system of coordinates ($,n> 5 ) shown In the inset In 
figure 3, Let (x,y,z) be the coordinates of a point on the propeller blade 
at r/R * 0.75 for a given azimuth angle. The column vector U represents 
the velocity components. Matrix A Is the t# rotation matrix and V 
represents the transformed vector. The transformation is given by 

V - AU (1) 

If B represents the ig rotation matrix, the final transformed vector Is W. 
This transformation Is given by 

W ® BV (2) 

The final transformed vector W Is related to U by 

W * BAU (3) 

Equation (1) may be written as 
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The ui, U 2 , and U 3 are the x, y, and z velocity components given by the 
method of reference 4 at the point (x,y,z). Equation (2) may be written as 
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Here the wi , W2 , and W3 are the velocity components that are shown in 
figure 3. 
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To verify the results of the present investigation, velocities at off- 
body points that were computed by the method of reference 4 were compared to 
those computed by two different methods, one of which was the transonic 
potential flow solution of Jameson (ref. 6 ) and the other was the modified 
small disturbance theory program (ref. 7). Since the Jameson method cannot 
treat wing-body configurations, a wing-alone case was computed using each of 
the three methods. The geometric characteristics of the wing alone is identi- 
cal to the wing of the PBW 4 N configuration from the 12% to the 100% semi- 
span stations. The wing area of the wing-alone configuration is approximately 
equal to the exposed wing area of the PBW 4 N configuration. The Jameson 

computer program was modified to print velocities at off-body mesh points in 

the sheared parabolic coordinate system. The mesh points were in a region in 
front of the wing, above and below the wing chord plane. The coordinates of 

the selected mesh points were used as inputs to the method of reference 4 

which has the capability of computing velocities at arbitrarily specified off- 
body points, so that a direct comparison of u .he velocities can be made. Like 
the Jameson program the method of i*eference 7 does not have the capability of 
computing velocities at arbitrary tff-body points. This computer program, 
however, was similarly modified to print velocities in a given region of the 
wing-alone flow field, Since the program has been designed to generate its 
own coordinate system, it was necessary to Interpolate between mesh points to 
obtain velocities at given "Jemeson mesh points," Shown in figures l 6 (a-c) 
are comparisons of the various velocity components. The Ax/c indicates the 
distance ahead of the wing leading edge. The coordinates have been normalized 
by the mean aerodynamic chord and the semispan of the PBW 4 N configuration, 

All three methods agree fairly well with one another with the exception of 
the method of reference 7 which predicts lower values of the w/V m than the 
other two methods. 

Shown in figure 17 are the velocities at off-body points corresponding 
to points in the propeller disc at r/R « 0.75 for the wing-alone cases 
computed by the three different methods. Note that, as in the previous com- 
parisons (fig. 16(c)), the overall level of the w/V ro component computed by 
the method of reference 7 is lower than the w/V^ levels computed by the 
methods of references 4 and 6 . The effect of the differences in w/V w on the 
blade angle-of-attack characteristics is shown in figure 18, The wing-alone 
velocities were adjusted for the effects of the body and nacelle using incre- 
ments computed by the method of reference 4. The blade angle-of-attack char- 
acteristics based on the velocities computed by the methods of references 4 
and 6 are shown to be in good agreement with each other while based on 

the results of reference 7 shows a different overall level. 


L 
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TABLE 1, 


AIRFOIL COORDINATES FOR WINGS Wj AND W 2 


0.00000 
. 00901 
,02153 
.03806 
.05904 
.08427 
,11350 
.14645 
.18280 
.22222 
.26430 
.30866 
.35486 
.40246 
,45099 
.50000 
.54901 
.59755 
. 64514 
.69134 
.73570 
.77770 
.81720 
.85355 
.88651 
.91574 
.94096 
.96194 
.99039 
1.00000 


Upper 

surface 

0.00000 

.02406 

.03579 

.04677 

.05650 

.06450 

.07045 

.07432 

.07638 

.07695 

.07635 

.07476 

.07231 

.06908 

.06520 

.06074 

.05579 

.05047 

.04490 

,03918 

.03345 

.02782 

.02243 

.01744 

.01297 

.00912 

,00597 

.00353 

,00067 

.00000 


Lower 

surface 

0.00000 

-.02406 

-.03579 

-.04677 

-.05650 

-.06450 

-.07045 

-.07432 

-.07638 

-.07695 

-.07635 

-.07476 

-.07231 

-.06908 

-.06520 

-.06074 

-.05579 

-.05047 

-.04490 

-.03918 

-.03345 

-.02782 

-.02243 

-.01744 

-.01297 

-.00912 

-.00597 

-.00353 

-.00067 

.00000 



TABLE 2.— AIRFOIL COORDINATES FOR TUNG W- 
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(a) Geometric characteristics. 
Figure 1.— Nacelle geometry. 
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(b) Pitch and yaw sign convention. 
Figure 1.— Concluded. 
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Figure 3.— Effect of the nacelle on the blade angle-of -attack 
characteristics of the propeller (P) ; 1 » 0°, 1. * 0°. 
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Figure 8,— Effect of nacelle Incidence of the blade angle-o£~attaek 
characteristics for the PBWiN configuration; 1,, « 0°. 
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